Mon. Not. R. Astron. Soc. 000, 1-12 (2012) Printed 10 October 2012 (MN KTpX style file v2.2) 



SN 2007bg: The Complex Circumstellar Environment Around One 
of the Most Radio-Luminous Broad-Lined Type Ic Supernovae 

P. Salas, 1 F. E. Bauer, 12 C. Stockdale 3 ' 4 J. L. Prieto, 5 * 

1 Pontificia Universidad Catolica de Chile, Departamento de Astronomia y Astrofi'sica, Casitla 306, Santiago 22, Chile 

2 Space Science Institute, 4750 Walnut Street, Suite 205, Boulder, CO 80301 

3 Marquette University, Physics Department, PO Box 1881, Milwaukee, WI 53201, USA 

4 The Homer L. Dodge Department of Physics & Astronomy, 440 W. Brooks ST, The University of Oklahoma, Norman, OK 73019, USA 
5 Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA 



10 October 2012 

ABSTRACT 

In this paper we present the results of the radio light curve and X-ray observations of 
^ <-h broad-lined Type Ic SN 2007bg. The light curve shows three distinct phases of spectral and 

Q_j temporal evolution, implying that the SNe shock likely encountered at least 3 different cir- 

cumstellar medium regimes. We interpret this as the progenitor of SN 2007bg having at least 
two distinct mass-loss episodes (i.e., phases 1 and 3) during its final stages of evolution, yield- 
ing a highly-stratified circumstellar medium. Modelling the phase 1 light curve as a freely- 
expanding, synchrotron-emitting shell, self-absorbed by its own radiating electrons, requires 
a progenitor mass-loss rate of M ps 1.9 x 10~ 6 (w M ,/1000 kms~ ) M Q yr _1 for the last 
t ~ 20(1^/1000 kms~ ) yr before explosion, and a total energy of the radio emitting ejecta 
of E ps 1 x 10 48 erg after 10 days from explosion. This places SN 2007bg among the most en- 
ergetic Type Ib/c events. We interpret the second phase as a sparser "gap" region between the 
two winds stages. Phase 3 shows a second absorption turn-on before rising to a peak luminos- 
ity 2.6 times higher than in phase 1. Assuming this luminosity jump is due to a circumstellar 
medium density enhancement from a faster previous mass-loss episode, we estimate that the 
phase 3 mass-loss rate could be as high as M < 4.3 x 10~ 4 (v„,/1000 kms -1 ) M Q yr -1 . The 
phase 3 wind would have transitioned directly into the phase 1 wind for a wind speed differ- 
ence of ps 2. In summary, the radio light curve provides robust evidence for dramatic global 
changes in at least some Ic-BL progenitors just prior (~ 10 — 1000 yr) to explosion. The ob- 
served luminosity of this SN is the highest observed for a non-gamma-ray-burst broad-lined 
Type Ic SN, reaching is. 46 GHz ~ 1 x 10 29 erg Hz -1 s _1 , ~ 567 days after explosion. 
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1 INTRODUCTION 

Core collapse supernovae (SNe) mark the spectacular death of mas- 
sive stars (Woosley & Weaver 1986). Among these SNe Types lb 
and Ic (collectively referred to as Ib/c due to their perceived similar- 
ities) are associated with the explosion of massive stars which were 
stripped of their H envelopes by strong stellar winds. Types lb show 
He features in their spectra, whereas Types Ic do not. Our current 
understanding of Type Ib/c SNe suggests that Wolf-Rayet (WR) 
stars are a primary progenitor path (Van Dyk, Li & Filippenko 
2003), although the relative numbers of WR stars fall well short of 
the observed Type Ib/c rates and thus a lower-mass binary progeni- 
tor path is most likely required to account for a significant fraction 

* Hubble, Carnegie-Princeton Fellow. 



of these events (e.g., Podsiadlowski et al. 1992; Ryder et al. 2004; 
Weiler et al. 201 1 ; Smith et al. 201 1). Archival imaging of pre-SNe 
locations, however, have failed to produce direct associations be- 
tween Type Ib/c SNe and their possible progenitors (e.g., Smartt 
2009). Thus we must rely both on theory and indirect information 
obtained after the explosion. Fortunately, X-ray and radio emis- 
sion are detectable by-products of the shock interaction between 
the high-velocity ejecta and the low-velocity progenitor wind in 
SNe (whereas the optical emission is generated by lower-velocity 
ejecta and radioactive decay). These observables crudely scale with 
the density of the circumstellar medium (CSM), and hence pro- 
vide constraints on the wind properties of the SN progenitor (e.g., 
Chevalier 1982a,b). Because massive stars are expected to have 
strong episodic mass-loss (e.g., Humphreys & Davidson 1994) and 
line driven winds (e.g., Vink, de Koter & Lamers 2000), radio mon- 
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itoring of Type Ib/c SNe provide interesting constraints on possible 
progenitor types. The past decade has witnessed increased attention 
on Type Ib/c SNe due to the association of Type Ic SN 1998bw with 
gamma-ray burst (GRB) 980425 (Galama et al. 1998), as these ob- 
jects may help elucidate how the central engines of both GRBs and 
their associated SNe work. Following this historic connection, sev- 
eral other Type Ic SNe have been linked to GRBs (e.g., Hjorth et al. 
2003; Kawabata et al. 2003; Matheson et al. 2003; Staneketal. 
2003; Frail et al. 2005; Modjaz et al. 2006; Woosley & Bloom 
2006; Bufano et al. 2012); all of these objects show broad emis- 
sion lines in their optical spectra, with velocities above 15000 km 
s - . To date so far 5 GRBs have been spectroscopically confirmed 
to have a SNe counterpart, while 14 show photometric evidence 
of a SNe counterpart (Hjorth & Bloom 2011), and hence a larger 
sample is needed in order to shed some light on the explosion 
mechanisms involved. Given the large number of SNe Type Ic com- 
pared to GRBs, one key question is whether these SNe could sim- 
ply be off-axis GRBs. Initial efforts concluded that the vast major- 
ity were not (e.g., Soderberg et al. 2006c). Recently new evidence 
was presented for the detection of at least one possible off-axis 
GRB with no direct observation of the gamma-ray emission: SN 
2009bb (Soderberg et al. 2010b); as well as the more controver- 
sial SN 2007gr (Paragi et al. 2010; Soderberg et al. 2010a; Xu et al. 
2011). 

Here we report on a particularly interesting broad-lined Type 
Ic, SN 2007bg, which was optically discovered on 2007 April 16. 15 
UT at a = ll h 49 m 26.18 s , 5 = +51°49'21.8" (J2000) at red- 
shift z = 0.0346 (Quimby, Rykoff & Yuan 2007) (d = 152 Mpc, 
H = 70 km s" 1 Mpc' 1 , Q M = 0.3, and fi A = 0.7 in a ACDM 
cosmology). A spectrum taken on 2007 April 18.3 with the 9.2 
m Hobby-Everly Telescope let Quimby et al. (2007) classify SN 
2007bg as a Type Ic broad-lined (Ic-BL) supernova. The classifica- 
tion of this SN as a Type Ic and the fact that it resides in a faint host 
encouraged Prieto, Stanek & Beacom (2008) to think about this ob- 
ject as a likely off-axis GRB. Prieto, Watson & Stanek (2009) also 
noted that SN 2007bg was one of the brightest radio SNe one year 
after explosion, making it an even better candidate for an off-axis 
GRB. However, based on the radio light curve of SN 2007bg known 
at the time, Soderberg (2009) contended that the observed radio 
emission could be explained by the presence of a density enhance- 
ment in the CSM, rather than a GRB like central engine driving the 
explosion. Using the archival radio data for SN 2007bg, we further 
address this issue below. 

The article outline is the following. In §2 we introduce the 
observations taken with the Very Large Array (VLA). Our results 
and modelling of the radio light curves are presented in §3. Finally 
in §4 we summarise our findings and provide a brief discussion on 
the future of radio observations of SNe. 



2 OBSERVATIONS 
2.1 Radio 

The VLA 1 radio observations were carried out between 2007 April 
19 and 2009 August 28 as part of the programs AS887, AS929 and 
AS983 (PI A. M. Soderberg). 

1 The Very Large Array and Very Long Baseline Array are operated by the 
National Radio Astronomy Observatory, a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universi- 
ties, Inc. 



The data were taken in standard continuum observing mode, 
with a bandwidth of 2 x 50 MHz. Data were collected at 1.43, 
4.86, 8.46, 15, 22.5 and 43.3 GHz. As primary flux density cali- 
brators 3C286 and 3C147 were used. We used these calibrators to 
scale our flux density measurements to the Perley-Buttler 2010 ab- 
solute flux scale. For secondary calibrators three different sources 
were used; SDSS Jl 14856.56+525425.2 in most epochs, as well as 
SDSS Jl 14644.20+535643.0 and SBS 1 150+497. The derived flux 
density for the secondary calibrators is shown in Table 1 . Tradition- 
ally one has been able to cross-check the derived secondary cali- 
brator fluxes with archival data, but in this case the VLA stopped 
regular monitoring of calibrators since 2007. This, in combination 
with the weak calibrator for the 15 and 22.5 GHz data 2 , makes this 
data less reliable than the rest at lower frequencies, even adopting 
larger error bars for it. 

During this period of time, the VLA incorporated the new 
EVLA correlator to some of the antennas, so some of the antennas 
did not have all of the different wavelength receivers. This greatly 
reduced the sensitivity of observations during some epochs, espe- 
cially at 1.43 and 15 GHz. We flagged, calibrated, and imaged 
the observations using the CASA' 5 (Common Astronomy Soft- 
ware Applications) software following standard procedures. We 
combined VLA and EVLA baselines, and to reduce the effects 
of the differences between VLA and EVLA antennas, we derived 
baseline-dependent calibrations on the primary flux calibrator. 

To measure the flux density of the source we fitted a Gaussian 
at the source position in our cleaned images, as for upper limits 
these were computed from the dirty map rms in the central region. 
Since the rms noise in each cleaned map only provides a lower 
limit on the total error we added a systematic uncertainty e„ due 
to possible inaccuracies of the VLA flux density calibration and 
deviations from an absolute flux density scale resulting in a final 
error for each flux density measurement of a 2 — (€„S V ) 2 + (J 2 ms , 
with e„ = 0.2,0.2,0.1,0.05,0.05,0.05 at 43.3, 22.5, 15, 8.46, 
4.86, and 1.43 GHz. The factors were chosen based on the quality 
of the secondary calibrator, and on the observed phase scatter after 
calibration of the data. The resulting flux density measurements of 
SN 2007bg are shown in Table 2 along with the image rms. 

2.2 X-Ray 

SN 2007bg was observed seven times with the Swift X-ray Tele- 
scope (XRT) and once with the Chandra X-ray Observatory, as 
listed in Table 3. Processed data were retrieved from the Swift and 
Chandra archives, and analysis was performed using the HEA- 
SOFT (v6.12) 4 and CIAO (v4.4.1) 5 software packages, along with 
the latest calibration files available at the time. In the case of Chan- 
dra, we reprocessed the data to apply the latest calibration modifi- 
cations, apply positional refinements, and correct for charge trans- 
fer inefficiency (CTI). All datasets were then screened for standard 
grade selection, exclusion of bad pixels and columns, and intervals 
of excessively high background (none was found). 

The seven Swift datasets were combined into three rough 
epochs. We adopted a 35"radius aperture (corresponding to ~85% 
encircled energy) as the source extraction region and modelled 
the local diffuse and scattered background using an annulus of 

2 http : / /www . via . nrao . edu /astro /calib/ manual 

3 http://casa.nrao.edu/ 

4 http : / /heasarc .gsfc.nasa. gov/lheasof t/ 

5 http://asc.harvard.edu/ciao/ 
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Table 1. VLA secondary calibrator measurements 3 



Obs. 


A.43 GHz 


/4.86 GHz 


f8. 46 GHz 


/l5 GHz 


/22.5 GHz 


Date 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


2UU / -Apr- 1 9 






0.4525 






ZUU /-Api-Zj 






A A Qrtt 

U.4.5UO 




A T7"21 
0.2/ jl 


ZUU / -Apr- 24 




0.4571 


0.433 






ZUU /-Apr-Zj 






0.453 




0.282 


inm An.- 
ZUU /-Api-Zo 




a a ceo 
0.4:D2 


0.4521 


0.3644 




inrn An.- id 
ZUU /-Api-JU 




a /i /C"7 
U.4o/ 


0.4583 






ZUU /-May -J 








0.355 


0.2752 


inm i j„,. c 

2UU/-May-5 




0.458 


0.4427 






ZUU / -May- 1 Z 




ft All 

U.41Z 


0.396 


U. JoU/ 


0.2733 


ZUU / -May- 1 / 




a a ^oo 
U.4joZ 


A A&11 
U.4oZ3 


A I^OA 

u. jjyu 


A 

U.Zoo 


ZUU /-May-Z / 


i\ a c n i b 

U.4jJj 


0.4507 


0.4173 


0.345 


0.2784 


20U7-Jun-1 1 




0.45 


0.464 


0.359 


0.2781 


taat Tuvi to 
ZUU/-Jun-ZZ 




0.4587 


0.4372 






taat t..i i 
ZUU /-Jul- j 








0.30 


0.29 


onm t..i i 
ZUU / -Jul- / 




0.4599 


0.4479 






ZUU/-JU1-Z4 


0.465 


0.453 


0.4462 


0.349 


0.266 


ZUU/ -Aug- lo 




a /i ^a ^ 


n a a 1 n 

U.441U 


U. J4o 


A T7TA 
U.Z /ZU 


iaat Can n 
ZUU /-iep- / 




a /i^/i 
U.4o4 


A AA~7^ 
U.44 / J 


A l^T 

U.JjZ 


A "2 A 
U.JU 


2UU/-Sep-23 


0.476 b 


0.4668 


0.4624 






ZUU /-Uct-Z-5 


0.459 b 


0.45 14 


0.4410 






ZUU/-JNOV-1 / 


n /i cto b 
U.4joZ 


0.4332 


A /I /I /I A 

U.444y 






ZUU /-Dec-Zj 


0.452 b 


0.4583 


A /I C 1 1 

0.43 1 J 






TAAQ Ton 1 

zuuo-Jan-j 










U.Zooj 


2008-Jan-28 




0.50 


0.475 




0.10 


2008-Feb-25 




0.4642 


0.4561 




0.304 


2008-Apr-19 


0.476 b 


0.6247 b 


0.587 b 






2008-May-7 






0.449 






2008-Jun-9 




1.197 c 


1.05 c 




1.050 c 


2008-Nov-3 


0.5160 b 


0.6480 b 


0.6069 b 




0.49 l b 


2008-Dec-30 






0.481 


0.351 




2009-Apr-6 


0.5100 b 


0.575 b 


0.5430 b 




0.3491 b 


2009-May-31 


0.503 b 


0.5763 b 


0.523 b 






2009-Aug-27 


0.530 b 


0.5411 b 


0.5162 b 







a Unless otherwise stated, the flux density measurements are 
for SDSS Jl 14856.56+525425.2. 

b Flux density measurements are for SDSS Jl 14644.20+535643.0. 
c Flux density measurements are for SBS 1150+497. 
On 2008 January 5 the source was observed at 43.3 GHz. 
The secondary calibrator flux during this epoch is 0.209 Jy. 



35-100". No statistically significant detection was found within 
our adopted aperture. For the Chandra dataset, we adopted a 
l'.'5 radius aperture (corresponding to ~95% encircled energy) 
as the source extraction region and modelled the local diffuse 
and scattered background using an annulus of 1.5-10". Although 
five (5) counts fell within our aperture (a ~ 2<r signal), no sta- 
tistically significant detection could be made. Upper limits in 
all cases were determined using the Bayesian technique from 
Kraft, Burrows & Nousek (1991). We used PIMMS v4.6 to deter- 
mine a count-to-flux conversion, adopting an absorbed APEC ther- 
mal plasma model with kT = 4 keV and Ah = 1.5 x 10 20 cm -2 
(Dickey & Lockman 1990). 



3 RESULTS 

The resulting light curves for SN 2007bg are shown in Fig. 1 . The 
evolution of the light curves can be divided in three phases. Phase 
1: this phase starts after the explosion and lasts until day 144. Here 
we observe the typical turn-on of radio supernovae (RSNe), first 



starting at shorter wavelengths and cascading to longer wavelengths 
with time as wavelength-dependent absorption becomes less im- 
portant. The peak luminosity reached during phase 1 is 4.1 x 10 28 
erg Hz -1 s _1 observed on day 55.9 at 8.46 GHz. This makes SN 
2007bg one of the most luminous Type Ib/c RSNe on comparable 
time-scales. The first detections at 4.86 and 8.46 deviate from the 
rest of the evolution, and it could be due to interstellar scintillation, 
a real deviation from pure SSA due to clumpy CSM, or decelera- 
tion; unfortunately there is not enough early-time data to confirm 
it robustly. On day 55.9 we also observe a rise in luminosity at 15 
and 22.5 GHz, this could be the effect of a small density enhance- 
ment, since that would produce the achromatic variation observed, 
and as the lower frequencies are still optically thick it would not 
be noticed at these frequencies. Phase 2: if the shock is expanding 
into a constant (or even smoothly-varying) stellar wind, then after 
the absorption turn-on one expects the source to fade in intensity 
following a power-law decline with time. However on day 144 the 
flux density consistently deviates from this at 8.46 and 4.86 GHz, 
showing a discontinuous drop in the observed flux density from 
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Table 2. VLA flux density measurements for SN 2007bg 


Obs. 


At 


A. 43 GHz /4.86 GHz 


/8.46 GHz 


/l5 GHz 


/22.5 GHz 


/43.3 GHz 


Array 


Date 


(Days) 3 


(/ijy) (/iJy) 


OtJy) 


(nJy) 


OJy) 


(MJy) 


Configuration 


2007-Apr-19 


2.9 




$264 








D 


2007-Apr-23 


7.0 




$348 




sg 1421 




D 


2007-Apr-24 


7.8 


s£390 


$387 








D 


2007-Apr-25 


8.8 




$417 




1293 ± 220 




D 


2007-Apr-26 


9.8 


^339 


$555 


$666 






D 


2007-Apr-30 


13.8 


sg291 


480 ± 102 








D 


2007-May-03 


16.8 






$978 


1032 ± 75 




D 


2007-May-05 


19.2 


^1692 


753 ± 107 








D 


2007-May-12 


26.1 


^1905 


804 ± 100 


1900 ± 300 


1062 ± 245 




D 


2007-May-17 


30.9 


^657 


728 ± 87 


1480 ± 289 


$1050 




D 


2007-May-27 


41.3 


171 487 ± 96 


1257 ± 53 


1250 ± 270 


$1020 




A 


2007-Jun-ll 


55.9 


^1191 


1490 ± 104 


2000 ± 308 


1246 ± 127 




A 


2007-Jun-22 


66.8 


709 ± 60 


1390 ± 45 








A 


2007-Jul-03 


77.8 






$843 


$978 




A 


2007-Jul-07 


81.8 


915 ± 62 


1325 ± 47 








A 


2007-Jul-24 


98.8 


^201 886 ± 80 


1131 ± 59 


$909 


$870 




A 


2007-Aug-18 


124.0 


1050 ± 90 


957 ± 64 


$957 


$693 




A 


2007-Sep-07 


144.0 


^336 1103 ±88 


621 ± 81 


^1 137 


$825 




A 


2007-Sep-23 


159.8 


^525 792 ± 77 


316 ± 57 








AnB 


2007-Oct-23 


189.9 


^585 685 ± 62 


379 ± 59 








AnB 


2007-Nov-17 


214.9 


^342 479 ± 65 


404 ± 62 








B 


2007-Dec-23 


250.9 


^387 685 ± 65 


783 ± 57 








B 


2008-Jan-03 


261.8 








1723 ± 77 




B 


2008-Jan-()5 


263.8 










$1662 


B 


2008-Jan-28 


286.8 


sJ594 


1669 ± 90 




$999 




B 


2008-Feb-25 


314.8 


1059 ± 70 


2097 ± 73 




2570 ± 160 




CnB 


2008-Apr-19 


368.8 


^822 1304 ±62 


2200 ± 59 








C 


2008-May-07 


386.8 




2852 ± 58 








C 


2008-Jun-09 


419.9 


2039 ± 88 


3344 ±81 




2023 ± 140 




DnC 


2008-Nov-03 


566.9 


$321 2812 ±60 


3897 ± 78 




1854 ± 128 




A 


2008-Dec-30 


623.8 




3891 ±78 


$2358 






A 


2009-Apr-06 


720.8 


746 ± 168 3390 ±136 


3842 ± 69 




1629 ± 115 




B 


2009-May-31 


775.8 


$1167 3460 ±80 


3641 ± 79 








CnB 


2009-Aug-27 


863.8 


$1503 3373 ± 84 


3408 ± 109 








C 



a Days since explosion assuming an explosion date of 2007 April 16.15 UT. 
All upper limits are 3cr. 



Table 3. X-Ray measurements for SN 2007bg 



Date 3 


Instrument 


Obsid 


Exposure 


0.5-8.0 


Fo. 5-8.0 




-^0.5-8.0 






(ks) 




Counts 


(erg s _1 cm" 


- 2 ) 


(erg s _1 ) 


2007-AP1--22 1 


Swift XRT 


00030920001 


9.4 


<4.6 


< 2.0 X 10" 


14 


< 5.5 X 10 40 


2007-Apr-30 2 


Chandra ACIS-S3 


7643 


34.0 


<11.2 


< 2.6 X 10" 


15 


< 7.2 x 10 39 


2008-Mar-05 3 


Swift XRT 


00030920002 


2.8 










2008-Mar-17 3 


Swift XRT 


00030920003 


4.2 


< 1 1-1 


< 4.4 X 10" 


14 


< 1.2 x 10 40 


2008-Apr-23 3 


Swift XRT 


00030920004 


2.9 










2009-Jun-09 4 


Swift XRT 


00030920005 


3.5 










2009-Jun-15 4 


Swift XRT 


00030920006 


3.1 


<6.9 


< 3.1 x 10" 


14 


< 8.6 x 10 40 


2009-Jun-17 4 


Swift XRT 


00030920007 


2.0 











a Different epochs were combined to increase SNR. The superscript on the date indicates which were combined. 



days ~ 144 — 250. Phase 3: the luminosity of SN 2007bg begins a 
second turn-on around day 261.8, becoming brighter than in phase 
1 and reaching a peak luminosity of 1 x 10 29 erg Hz -1 s" 1 at 8.46 
GHz on day 566.9. This late time behaviour is comparable to that 
of other Type Ib/c SNe like SN 2004cc, SN 2004dk and SN2004gq 
(Wellons, Soderberg & Chevalier 2012), albeit much stronger and 
more well-defined than these other objects. There are also some 
Type lib's SNe which show late time modulations in their radio 



light curves. One of these is SN 2001ig (Ryder et al. 2004) in 
which case the late time modulations could be explained by a bi- 
nary companion of the progenitor. Another example is SN 2003bg 
(Soderberg et al. 2006a), where the late time variations can not 
rule out a binary companion or different mass-loss episodes of the 
progenitor. In both cases the modulations in the light curve result 
from modulations in the CSM density. A comparison between SN 
2007bg, other SNe and GRBs is shown in Fig. 2. During phase 1 
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Figure 1. Light curves for SN 2007bg. Phase 1 is shown as the first white region, phase 2 is the shaded region and phase 3 the last white region. Different 
frequencies are represented by the following colours: 1.43 GHz, 4.86 GHz, 8.46 GHz, 15 GHz, 22.5 GHz and 43.3 GHz. We see a slow turn-on at 8.46 
GHz during phase 1. The lines during phase 1 represents the model fits. The thin line assumes a density profile r oc t 094 and B oc t~ \ while the thick line 
requires r oc t 0,75 and B oc t~ 0,75 . During this phase we modelled the emission as arising from a spherical expanding shell self that is absorbed by its own 
radiating electrons, as explained in §3.2. During phase 3 the line model represents our best fit model including the presence of SSA and a clumpy CSM. This 
is described in §3.4. Filled circles denote detections with lcr error bars (which are smaller than symbols in some cases), while inverted triangles denote 3<r 
error limits. 



the light curve of SN 2007bg resembles that of other SNe. It shows 
the characteristic turn-on, but with a higher luminosity than other 
normal Ib/c SNe and comparable to GRB-SNe. Then during phase 
3 we observe a strong rise in its luminosity, which is unprecedented 
for SNe, as it reaches such high luminosities, 2.6 times larger than 
during phase 1 . 

The spectral index, /„ oc v a , evolution of the source is shown 
in Fig. 3. We observe a trend towards a negative value between the 
best sampled frequencies 4.86 GHz and 8.46 GHz during phase 1, 
but it could be that during this phase the source is still optically 
thick at 4.86 GHz. If we assume that at this time the emission has 
reached its peak luminosity at 4.86 GHz, then the late observed 
values of the spectral index point to a value of a ~ —1. Even 
the deviant points at 15 and 22.5 GHz on day 55.9 are consistent 
with this spectral index value. During phase 2 the source spectrum 
begins to invert. This behaviour continues during the first days of 
phase 3, indicating that some absorption mechanism becomes im- 
portant again in phase 3. From the earliest spectral indices avail- 



able we observe that the values are a ~ 2, which is lower than 
the typical synchrotron self absorption (SSA) value of 2.5 and the 
spectrum is shallower than the exponential cut-off expected from 
free-free absorption (FFA). This suggests may be observing a va- 
riety of emission regions from SN 2007bg; some being absorbed 
more highly than others. Since we can not resolve the source, we 
get the combined flux resulting in a shallower optically thick part 
of the spectrum. 

Given that we do not have strong constraints on how the intrin- 
sic spectral index varies with time, we adopt a fixed spectral slope 
for each phase of a ~ — 1. 

3.1 Early Physical Constraints from Phase 1 

Our first data point at 8.46 GHz allows us to evaluate the bright- 
ness temperature of the source at 13.8 days after explosion. The 
brightness temperature of the source should not exceed 3 x 10 11 
K if it is to remain under inverse Compton catastrophe (ICC), 
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Figure 2. Radio light curves for different Type lb SNe (shades of blue with different symbols) Type Ic's (shades of pink with different symbols), a Type lib 
SN (green) and GRBs (black with different symbols). SN 2007bg (red stars) is shown along with other objects. We observe that the radio luminosity from SN 
20()7bg is lower than that of GRBs on comparable time-scales and it is larger than in regular SNe at the time of the first peak during phase 1. Phase 3 starts 
at a comparable time-scale compared to other SNe, but the peak luminosity reached is larger by orders of magnitude. The SNe shown in this plot are Type Ic 
2002ap (circles; Berger, Kulkarni & Chevalier 2002), Type lib 2003bg (circles; Soderberg et al. 2006a), Ic 2003L (squares; Soderberg et al. 2005), Ic 2004cc, 
lb 2004dk, lb 2()04gq (diamonds; circles; squares; Wellons et al. 2012), Ic 19941 (pentagons; Weiler et al. 201 1), Ic 2007gr (triangles; Soderberg et al. 2010a), 
lb 2007uy, Ib/X-ray flash 2008D (diamonds; triangles; van der Horst et al. 2011) and Type Ic broad-lined 2009bb (hexagons; Soderberg et al. 2010b). Also, 
SNe associated with a GRBs Type Ic 1998bw (circles; Kulkarni et al. 1998), Ic 2003dh (squares; Hjorth et al. 2003; Kawabata et al. 2003; Matheson et al. 
2003; Frail et al. 2005), Ic 20031w (diamonds; Soderberg et al. 2004) and Ic 2006aj (Xs; Soderberg et al. 2006b). 



where the cooling time due to inverse Compton scattering is shorter 
than the synchrotron cooling time (Readhead 1994) (see also; 
Kellermann & Pauliny-Toth 1969). Using the measured flux den- 
sity of 753 /iJy at 8.46 GHz, and if we assume that the optical ex- 
pansion speed of ~ 17500 km s~ , measured from the optical Ca II 
triplet, at ~ 5 days after explosion (Young et al. 2010) also pertains 
the radio emitting region, we obtain T « 3.6 x 10 12 K. This value 
is above the ICC limit, so one of our assumptions must be wrong. 
As discussed in the introduction, typically the radio ejecta from 
SNe have higher velocities than that of the photosphere, which pro- 
duces the optical emission. Based on this we adopt a radius for the 
source such that the brightness temperature remains below the ICC 
limit constraining the source size to be (7 ± 1) x 10 cm. The re- 
ported errors come from the flux density measurement uncertainty. 
At 13.8 days after explosion this implies a bulk expansion speed of 
/3F w (0.2 ± 0.04)c, where ft = v/c and T is the Lorentz factor. 
This speed v ~ (0.19 ± 0.03)c is ~ 3 times larger than that of 
the photosphere. We stress that this only a lower limit on the actual 
velocity at which the radio emitting region expands, as the radio 
ejecta could be expanding at a higher velocity and be well below 



the ICC brightness temperature. The resulting brightness tempera- 
ture is shown in Fig. 4 where we assumed that the radio emitting 
region does not decelerate i.e., r oc t. 

Another way to estimate the expansion speed of the radio 
emitting region can be obtained by placing SN 2007bg in a lu- 
minosity peak time plot, which as discussed by Chevalier (1998) 
provides information about the average expansion speed of the ra- 
dio emitting region. The plot is shown in Fig. 5 and we note that SN 
2007bg has an average expansion speed of ~ 0.2c typical of Type 
Ib/c RSNe and it shows one of the highest luminosities among Type 
Ib/c's and comparable to those of GRBs. It is important to note 
that other absorption mechanisms, such as FFA, or deviations from 
equipartition, would increase the expansion velocity (and hence 
source energy, derived later). 



3.2 Constraints from SSA Evolution during Phase 1 

Synchrotron emission from SNe originates from the shock region 
between the SNe ejecta and the pre-existing CSM. This shock am- 
plifies the existing magnetic fields via turbulence. The electrons 
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Figure 3. Spectral indices between adjacent bands. Over plotted are the predicted values from the light curve fitting as described in §3.2 and §3.4. The thin and 
thick lines represent the same models as in Fig. 1 . Top left panel shows the evolution of the spectral index between 20 cm and 6 cm. Top right panel shows the 
spectral index from 6 cm and 3.6 cm, the best sampled set of observations. Bottom left spectral index between 3.6 cm and 2 cm. Bottom right spectral index 
between 2 cm and 1 . 2 cm. We note that near the end of phase 1 a large negative spectral index is observed, this is mainly because the light curve at 1 . 2 cm 
deviates from the evolution observed at other wavelengths. Filled circles denote detections with lcr error bars (which are smaller than symbols in some cases), 
while inverted triangles denote 3<r upper limits and triangles 3<r lower limits. 



in the shock region then interact with this magnetic field produc- 
ing the observed synchrotron emission. This emission is subject to 
self-absorption, which is caused by SSA (Chevalier 1998) in the 
case of low density environments. To describe the time evolution 
of the flux density during phase 1 we adopt the parametrization 
of Soderberg et al. (2005) in which the radio emission originates 
from a spherical thin shell of width r/r/ expanding at a velocity v 
in which the intrinsic synchrotron emission is absorbed by SSA. 
The parameter r\ is the number of shells that fit into the size of the 
source with a value of 77 ^ 10 (Frail, Waxman & Kulkarni 2000). 
We adopt the assumptions of the standard model of Chevalier 
(1982b) where the hydrodynamic evolution of the ejecta is self- 
similar across the shock discontinuity. With this assumption the ra- 
dius of the source evolves as r oc t ar where the parameter a r is 
related to the density profile of the outer SN ejecta p oc r~ n and 
the density profile of the radiating electrons in the shock region 
n c oc r~ s such that a r = (n — 3)/(n — s). The self similar model 
also requires that s < 3, n > 5 and a r in the range 0.75 — 1. The 
model also assumes that the magnetic energy density and the rela- 
tivistic electron energy density scale as the total post-shock energy 
density. This implies that the fraction of energy contained in the 
relativistic electrons e c and magnetic field ee remains fixed. This 



in turn relates the magnetic field evolution B oc t" B , to the ra- 
dius of the source and the density profile of the radiating electrons 
a.B = (2 — s)a r /2 — l. This means that for a freely expanding wind 
s = 2, the magnetic field evolves as B oc t . For our analysis we 
adopt a source in equipartition with e c = eb = 0.1. Equipartition 
of energy minimizes the total energy. 

Using this model and requiring the time evolution of the 
magnetic field to go as B oc t" 1 we find that the radius of 
the source evolves as r oc t ' 94 . We will refer to this as model 
1. This is consistent with a shock expanding freely through the 
CSM. The normalization constants for the flux density and SSA 
are C f « 8.4 x 10~ 53 g s _1 and CV ~ 8.5 x 10 36 s 3 ' 5 , giv- 
ing a reduced chi square value of \r = 7.3 for 22 degrees of 
freedom. These values require the size of the source to go as 



cm, which in turn means that the ex- 

-0.06, 



r « 6.5 x 10 lj (i/i ) 

pansion velocity is v ~ 0.25(£/io) -u ' UD c. The magnetic field goes 
asBfs 2.6^/tor 1 G and the energy is E w 9.6 x 10 47 erg. For 
this we assumed a synchrotron characteristic frequency of ^ m = 1 
GHz. Using these values we can also estimate the mass-loss rate 
of the progenitor. Following Soderberg et al. (2006a) we estimate a 
mass-loss rate of M « 1.5 x 10" 6 (n„/1000 kms" 1 ) M Q yr _1 . 
This mass-loss rate lies at the low end of those observed for WR 
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stars in low metallicity environments (Crowther et al. 2002). The 
derived mass-loss rates from this model depend on the fraction of 
energy contained in relativistic electrons and magnetic fields, devi- 
ations from equipartition would require higher mass-loss rates and 
energy requirements. 

The early data hint at a shallower slope of the light curve, so 
we also fit a second model, model 2, where we allow the magnetic 
field evolution to change as B oc i -0 ' 75 and the source size as 
r oc t ' 75 . In this model the reduced chi square value is \r = 4.4 
for 21 degrees of freedom. In this case the normalization constants 
areC/ sa 1.7x 10~ 52 g s" 1 and CV f»4.6xl0 36 s 3 ' 5 . This implies 
a source size evolution of r w 8.7 x 10 15 (t/to) ' 74 cm, and an 
expansion velocity of v ~ 0.3(t/to)~°' 25 c. Given this expansion 
speed the source would remain below the ICC limit. The magnetic 
field goes as B ~ 2.1(£/to) -0 ' 75 G. The source energy with these 
values is E ~ 1.5 x 10 48 erg, and the mass-loss rate is M fti 
1.9 x 10 _6 (w ro /1000 kms" 1 )M yr" 1 .Here again we have used 
the same parameters for u m and 77. Given this time evolution for 
the magnetic field and source size, the density radial profile goes 
as p oc r~ 8 and the number density of emitting electrons as n c oc 
r -1 ' 3 . We adopt the result from model 2, as it provides a better fit 
to the data, although we cannot truly rule out model 1 due to the 
lack of early detections. 

Given this mass-loss rate, the circumstellar density for a wind 
stratified medium A*, defined as M /A-kv w — 5 x lO 11 ^* g cm -1 
(Li & Chevalier 1999), is A, w 0.2. Under these conditions a rel- 



ativistic jet from a GRB would become non-relativistic, and hence 
observable on a time-scale of £nr ~ 0.3(.Ek/10 51 erg)Aj 1 . With 
the kinetic energy of the ejecta of Ek ~ 4 x 10 51 erg, de- 
rived from the optical light curve (Young et al. 2010), this time- 
scale is £nr ~ 6 yr. The lack of clear afterglow signatures (e.g., 
Li & Chevalier 1999; Waxman 2004; Soderberg et al. 2006c) in 
our observations thus provide strong evidence against an expand- 
ing GRB jet being the source of the observed luminosity for SN 
2007bg. 



3.3 Stratified CSM: Phase 2 

As can be seen in Fig. 1, at around day 144 the flux density at 8.46 
and 4.86 GHz suddenly decreases creating a break in the power 
law flux decay of the optically thin radio emission. At 8.46 GHz 
the observed flux is a 44% less than the predicted one at the be- 
ginning of phase 2 on day 159.8. Based on the derived size of the 
radio emitting region, from the SSA modelling, this drop occurs at 
a distance of w 6.7 x 10 16 cm. The most likely cause for this de- 
crease in the radio emission is a drop in the CSM density. Here we 
consider two scenarios that would produce a change in CSM den- 
sity such as a binary companion or two different wind components 
from the progenitor star. In the first scenario, there should be evi- 
dence of this interaction from the optical spectra as Balmer series 
recombination lines. The late time optical spectra for SN 2007bg 
was too faint to be detected. The early optical spectra evidence 
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Figure 5. Source peak radio luminosity versus time of peak. The resulting position in this plot provides an estimate of the average expansion speed of the 
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points to a broad-lined SNe, most likely associated with SNe like 
1998bw (Young et al. 2010), and the latest optical spectra of SN 
2007bg where it is detected does not show evidence for a transi- 
tion from a type Ib/c to a type II. For this kind of SNe the most 
likely progenitor is a single WR. Also if there was a binary com- 
panion we would expect to see coherence between the emission 
before and after the drop/increase (depending on whether the CSM 
density dropped/increased at the region of interest). In the case of 
SN 2007bg, we see a > 1 dex jump in the normalizations of the 
unabsorbed phase 1 and phase 3 light curves, as well as a strong 
enhancement in absorption, making a binary companion CSM sce- 
nario highly unlikely. The second scenario considered is that the 
drop in flux density is due to a change in the progenitor wind. Based 
on the distance of the drop and with our assumed wind speed of 
~ 1000 km s _1 the change occurred ~ 45 yr before the explo- 
sion and lasted for « 23 yr. Since L v oc (M/«„) (a - 8+12a ' )/4 
(Chevalier 1982b) a flux density decrease of 44% implies that the 
mass-loss rate/wind speed dropped by a factor of ~ 5. For this es- 
timate we have used a r — 0.94. Our modelling of the light curves 
during phase 1 suggests that the density profile of the CSM is shal- 



lower (i.e., smaller a T ), but as discussed bellow, this relation over- 
estimates the mass-loss rate in such cases. 



3.4 Second Rise: Phase 3 

Phase 3 begins around day 261.8 with an increase in the flux den- 
sity at 22.5, 8.46 and 4.86 GHz. As discussed in §3 the luminos- 
ity of the source during this phase is larger than in phase 1 in all 
of the bands where we detect it. This flux density rise shows the 
characteristics of an absorption turn-on. At 22.5 GHz we observe 
what appears to be the peak of the turn-on on day 314.8, and at 
longer wavelengths the radio spectrum begins to turn-on near the 
end of observations. This turn-on is accompanied by an increase in 
the spectral index suggesting that the source has become optically 
thick again. The luminosity increase implies a sharp CSM density 
enhancement, presumably due to evolution of the progenitor wind, 
while the observed re-absorption requires either SSA or FFA. Fol- 
lowing Fransson, Lundqvist & Chevalier (1996) we can estimate 
the mass-loss rate from the progenitor assuming that the absorp- 
tion is purely FFA. Adopting the time of peak at 4.86 GHz as the 
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Table 4. SN 2007bg best-fit values for phase 3. 



Parameter 


Adopted value 


(X 


— l a 


p 


-l b 




l c 


Parameter 


Best-fit value 


Ki 


5 X 10 6 


K 3 


2.3 X 10 5 


6' 


-3.53 


K 5 


4.5 x 10 8 


5" 


-2.96 



a Adopted from observed values of a. 
b Assumed because of the lack of data. 
c Computed using a r = — (a — — 3)/3. 



time when the free-free optical depth is unity, and taking the expan- 
sion velocity of the shock from phase 1 then the mass-loss rate is 
M w 10 _1 (w l „/1000 kms -1 ) M Q yr" 1 . This rate is suspiciously 
high and inconsistent with our estimate (later in this section) of the 
mass-loss rate based on the unabsorbed portion of the spectrum. 
This high rate favours SSA. The absorption slope is relatively flat, 
and thus pure SSA with a similar model as fit to phase 1 does a poor 
job of fitting the turn-on. One possibility is differential SSA arising 
from small dense clumps in the CSM, which are preferentially af- 
fected by SSA while sparser regions are less affected; this naturally 
flattens the optically thick slope. The existence of clumpy winds in 
WR stars has been confirmed by observations (Moffat et al. 1988), 
and thus provides a natural way of explaining the observed radio 
emission. Alternatively, the density jump would naturally produce 
strong X-ray emission, which could ionize CSM material in front of 
the shock and provide the necessary free electrons for FFA. How- 
ever, without constraints on the X-ray luminosity and phase 3 CSM 
density distribution for SN2007bg, it is difficult to assess the feasi- 
bility of FFA. Given the poor fit of the pure SSA model (Xr = 9), 
we simply adopt the parametrization of Weiler et al. (2002) includ- 
ing only a clumpy CSM with the presence of internal SSA. In this 
case the flux density is normalised by K\, the clumpy CSM ab- 
sorption by Ks and SSA by Ks. We fit the above model to the light 
curve during phase 3. The best-fit values are given in Table 4. The 
resulting reduced chi square of this fit is \ r — 4.4 for 18 degrees 
of freedom. 

As shown in Fig. 1, this model provides a reasonable fit to the 
light curves at all frequencies and also to the spectral index, Fig. 3. 
The deviant point at 1.43 GHz could be easily affected by the host 
galaxy luminosity. The data did not allow us to constrain the decel- 
eration rate, so we adopted a constant expansion by fixing ft = — 1, 
which in turn requires that a r = m = 1 given the spectral index. 
Using our best-fit values we can estimate the mass-loss rate from 
the progenitor during this period following equations (11) and (13) 
of Weiler et al. (2002). The effective optical depth at the time of our 
latest observation, ~ 863 days, at 4.86 GHz is (r^ 5 ) « 2.2 x 10~ 3 
and m « 1. If we extrapolate the velocity derived in §3.2 to phase 
3, Vi « 0.13c at ti ~ 300 days, then the derived mass-loss rate is 
M « 4.3 x 10~ 4 (« ro /1000 kms -1 ) M yr~\ where an electron 
temperature of 10 4 K was used. This mass-loss rate is higher than 
the one observed during phase 1, which is a natural explanation 
for the high luminosity observed during phase 3. However, such a 
large mass-loss rate has not yet been directly observed, and it is 
at odds with the metallicity mass-loss rate relation derived from 



the Milky Way and Magellanic Clouds WR samples (see for e.g., 
Crowther 2006). This relation predicts orders of magnitude lower 
M in such a low metallicity host. On the other hand, this large 
mass-loss rate is also roughly consistent with that of the saturation 
limit for line-driven winds, M ~ 10~ 4 Mq yr _1 , and is of the 
same order of magnitude as values found in other Type Ib/c SNe 
(e.g., Wellons et al. 2012). We should caution however that the dra- 
matic flux increase between phases 1 and 3 implies a strong density 
contrast that is inconsistent with one of the primary assumptions of 
the Chevalier (1982b) self-similar model (namely a CSM density 
profile with a power-law slope). This density discontinuity could 
naturally lead to an overestimate of the mass-loss rate; to properly 
account for this type of case, one should generally model the hy- 
drodynamics in detail; unfortunately, such simulations are not well- 
constrained by radio data alone due to ambiguities in how the B- 
field and electron density distribution might change across the den- 
sity discontinuity and are beyond the scope of this work. Thus this 
phase 3 mass-loss rate should be regarded with caution and viewed 
as a likely upper limit. The strong mass-loss variations on such 
short time-scales seen here imply that the progenitor of SN 2007bg 
may have undergone unstable mass-loss just prior to explosion or 
a dramatic change in wind properties. Such rapid variations have 
been seen in a handful of objects before. For instance, in the case of 
Type lb SN 2006jc (Pastorello et al. 2007), a bright optical transient 
was observed two years prior to explosion. This was interpreted as 
a luminous blue variable-like outburst event from the progenitor, 
which the later SN explosion surely interacted with. Then there are 
other possibilities like the case of Type Iln SN 1996cr (Bauer et al. 
2008). In this event, the radio light curve also exhibits a strong late- 
time rise in emission associated with a density enhancement likely 
arising from the interaction of the SN blast-wave with a wind- 
blown bubble (Dwarkadas, Dewey & Bauer 2010). The progenitor 
of SN 1996cr is proposed to be either a WR or blue supergiant star. 
We still know relatively little about why such global changes in the 
progenitors occur, but it seems clear that a subset of SNe progeni- 
tors have them. 

3.5 X-Ray Constraints 

X-ray emission from supernovae can arise from three mechanisms; 
non-thermal emission from synchrotron radiating electrons, ther- 
mal bremsstrahlung emission from the material in the circumstellar 
and reverse shocks, and inverse Compton emission. Since X-ray 
emission from SN 2007bg could not be detected it is difficult to 
assess the feasibility of each of these mechanisms, so we will com- 
pare the derived mass-loss rates from the previous subsections with 
that predicted by the thermal bremsstrahlung process. 

If the X-ray emission comes from synchrotron radiating elec- 
trons then the X-ray luminosity can be obtained by extrapolating 
the radio synchrotron spectrum up to higher frequencies. To do this 
we must consider that the synchrotron spectrum changes its slope 
at the cooling frequency v c oc t~ 2 73 -3 bya factor Aa ~ —0.5 
(Soderberg et al. 2005). Using the derived magnetic field B from 
our SSA analysis of the radio light curve §3.2 we find that v c ~ 323 
GHz. Extrapolation of the luminosity as L„ oc v~ ' a from the cool- 
ing break on gives a X-ray luminosity 5 orders of magnitude bellow 
our upper limit ~ 30 days after explosion (epoch 2). Extrapolation 
of the cooling frequency up to epochs 3 and 4 of the X-ray data 
yields X-ray fluxes ~ 5 orders of magnitude bellow the upper lim- 
its. However, as mentioned in previous sections, the CSM disconti- 
nuity between phases 1 and 3 could cause a break in the magnetic 
field evolution. 
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X-ray emission from thermal bremsstrahlung comes from the 
interaction of the shocked electrons with the circumstellar material. 
This interaction heats the electrons in the medium which then cool 
via free-free emission. In this process the X-ray luminosity is re- 
lated to the density of the emitting material as L oc M 2 v^ u 2 t~ 1 
(Chevalier & Fransson 2003; Sutariaetal. 2003). We can obtain 
an upper limit on the mass-loss rate from the X-ray luminosity 
upper limits using equation (30) of Chevalier & Fransson (2006). 
At t ~ 30 days the X-ray luminosity is L . 5 _ g .o < 7.2 x 10 39 
erg s~ . We use the density profile of the outer SNe ejecta de- 
rived from our SSA analysis, p oc r" 8 . Given this luminosity and 
CSM properties thermal bremsstrahlung requires a mass-loss rate 
of M < 4.6 x 10~ 4 («„/1000 kms" 1 ) Moyr" 1 . At epochs 3 
and 4 of the X-ray data we obtain mass-loss rates of M < 3.6 x 
10~ 3 (« u ,/1000 kms-^andM < 1.4 x 10" 2 (^/1000 kms" 1 ) 
respectively. These mass-loss rates are consistent with the one re- 
quired by SSA during phase 1 and the light curve during phase 3. 

Inverse Compton scattering of the optical photons by the 
relativistic synchrotron emitting electrons can also be a cause 
of X-ray emission in SNe. It has been shown that this pro- 
cess played an important role in the case of SN 2002ap 
(Bjornsson & Fransson 2004), SN 2003L (Soderberg et al. 2005) 
and 2003bg (Soderberg et al. 2006a). Using equation (15) of 
Bjornsson & Fransson (2004) we obtain a ratio of U^/Ub ~ 12 
~ 30 days after explosion. For this we have used a bolometric lu- 
minosity of Lboi ~ 10 41 ' 5 erg s" 1 based on the pseudo-bolometric 
light curves for SN 2007bg of Young et al. (2010). Given the ratio 
U ph /UB and the radio flux of vF v k2x 10" 16 erg s" 1 cm" 2 the 
derived X-ray luminosity of~2xl0" 17 erg s" 1 cm -2 is in agree- 
ment with our upper limits. For later epochs we can not estimate the 
bolometric luminosity of SN 2007bg. 



4 SUMMARY AND CONCLUSIONS 

In this article we presented the radio light curves and X-ray ob- 
servations of Type Ic-BL SN 2007bg. The radio emission is char- 
acterised by three different phases. Phase 1 being characterised 
by a SSA turn-on. Based on the brightness temperature of SN 
2007bg the expansion speed of the radio emitting region is v — 
0.19 ± 0.02c. An expansion speed similar to what has been found 
in the small sample of radio-emitting Type Ib/c SNe. 

The derived mass-loss rate from the progenitor during phase 
1 is M « 1.9 x 10" 6 (i> m /1000 kms" 1 ) M Q yr"\ a low mass- 
loss rate among the observed in WR stars. During phase 2 we ob- 
serve a drop in the flux density while the spectral index seems to 
follow a smooth evolution. This drop is presumably due to the ra- 
dio ejecta entering a region of lower CSM density, probably aris- 
ing because of a change in the mass-loss rate of the progenitor 
or a variable wind speed. Then during phase 3 we observe a rise 
in the radio flux density along with a re-absorption of the radio 
emission. The shallow observed slope of the re-absorption could 
be explained by small clumps of material in the progenitor wind 
which would become more strongly self-absorbed than the sur- 
rounding CSM when shocked and produce a shallower optically- 
thick spectrum. From our modelling of the radio light curves dur- 
ing phase 3 we put an upper limit on the mass-loss rate of M < 
4.3 x 10" 4 (u™/1000 kms" 1 ) M yr" 1 . This second wind com- 
ponent could arise from a different stellar evolution phase of the 
progenitor prior to explosion, possibly a LBV phase. Or it could be 
the effect of stellar rotation on the stellar wind properties. If the ra- 
dio flux density from SN 2007bg declines in time like /„ oc t . It 



should still be detectable, with a flux density at 22.5 GHz of ~ 500 
/j,]y, and higher at lower frequencies. If the progenitor star went 
through more mass-loss episodes, then the observed flux density 
would deviate from this value. Recent observations taken with the 
Karl G. Jansky Very Large Array (JVLA) should help to establish 
this. Very few broad-lined type Ic SNe have radio detections, and 
none of them show signs of a complex CSM like that of SN 2007bg. 

These different wind components, along with the derived 
mass-loss rates will eventually help to constrain the progenitors of 
Type Ib/c SNe, from their inferred pre-SNe properties. Currently 
two main evolutionary paths has been proposed to account for the 
rate of these explosions: a single WR origin, and a binary system. 
However, the single WR scenario itself still needs significant clar- 
ification. For instance, what effects do evolution (mass-loss, rota- 
tion, etc.) and metallicity have in shaping the CSM around massive 
type Ib/c progenitors and ultimately their explosions. And more 
generally, how does this fold into the connection between type IPs, 
type Ib/c, and GRBs. 

Compounded by their relative rarity, a key impediment to im- 
proving our knowledge of broad-lined Ic SNe has been the lack 
of high-quality, multi-frequency observations from early through 
late epochs. Sensitive high-frequency radio observations of Ic SNe 
within 1-2 days of explosion are needed to efficiently identify the 
best objects for follow-up and allow further refinement of the phys- 
ical properties surrounding these unique explosions. Both ALMA 
and the newly retooled JVLA will play critical roles here, pro- 
viding rigorous probes of the synchrotron-emitting shock (to con- 
strain evolution of the shock velocity, any potential beaming, etc.) 
and ultimately new insights into the properties of their CSM (such 
as density, variability, dumpiness, and perhaps even dust content 
for these massive systems). Sensitive, long-term X-ray and opti- 
cal/NIR spectroscopic follow-up are also needed in order to break 
potential degeneracies and provide consistency checks against the 
shock-CSM interaction models. If such complete datasets can be 
assembled in the next several years, they may afford us a leap in 
our physical understanding of SNe and GRBs. 
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